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DYNAMIC  STALL  OF  AN  AIRFOIL 
WITH  LEADING  EDGE  BUBBLE  SEPARATION 
INVOLVING  TIME  DEPENDENT  RE-ATTACHMENT 

by 
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ABSTRACT 


The  dynamic  stall  of  an  airfoil  with  leading  edge 
bubble  separation  is  analyzed.  The  stall  flutter  of 
turbomachine  blading  often  involves  periodic  growth  and 
collapse  of  such  a bubble.  The  mathematical  model  rep- 
resenting the  physical  problem  is  presented.  A flat 
plate  undergoing  harmonic  oscillations  with  time  depen- 
dent point  of  re-attachment  is  studied  for  the  pertur- 
bed aerodynamic  reactions  and  applications  to  the  stall 
flutter  problem. 


♦Graduate  Research  Assistant  and  **Professor  of  Mechanical  Engineering 
Stevens  Institute  of  Technology 


Nomenclature 
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a = acceleration 
b = blade  semichord 
CL  = lift  coefficient 

CM  = moment  coefficient 

’f  - Fourier  transform 
f ,g=  auxiliary  functions 
h = translational  displacement 

i = imaginary  unit  used  in  representation  of  the  complex 
variables  Z and  £ 

j = imaginary  unit  used  in  representation  of  the  complex 
form  of  simple  harmonic  function 

k = reduced  frequency 

L = perturbation  lift 

M = perturbation  moment 

p = perturbation  pressure 

Rq  + j JQ  = complex  function  defined  by  Eq.  (14) 

+ j = complex  function  defined  b-  (15) 

S(t)  = point  of  reattachment 

S^/  S2  = upper  and  lower  limits  of  reattachment 
t = time 

u = perturbation  velocity  in  x-direction 
v = perturbation  velocity  in  y-direction 
V = freestream  velocity 

W = Q+iV,  complex  acceleration  potential 

x,  y = cartesian  coordinates  in  physical  plane 

Z = x+iy  point  in  physical  plane 

a = rotational  displacement  \ 

S = C+in  point  in  transformed  plane 
p = dummy  variable 

£ = coordinates  in  transformed  plane 

p = density  of  the  fluid 

a = a variable 


J 


$ *=  acceleration  potential 
X = phase  difference 
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Introduction 

An  important  type  of  phenomenon  classified  under  the  general 
heading  of  aeroelasticity  of  lifting  surfaces  is  stall  flutter. 

This  is  a type  of  dynamic  instability  that  occurs  when  the  flow 
separates  around  an  airfoil  through  the  whole  or  part  of  each  cycle 
of  its  vibratory  motion.  This  nonclassical  type  of  flutter  may 
involve  periodic  breakaway  and  reattachment  of  the  flow  and  various 
types  of  time  lag  effects  between  the  motion  and  the  reactions. 

The  stall  flutter  of  helicopter  rotors,  aircraft  wings,  aircraft 
engine  compressors  at  ground  start  up  and  high  speed  flight  are 
examples  giving  rise  to  non-steady  flow  about  an  airfoil  and  a 
possible  condition  of  destructive  behavior.  Thus  the  prediction  of 
aerodynamic  reactions  in  such  flow  situations  is  of  considerable 
importance . 


* 
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In  many  instances  the  flow  separation  from  the  suction  surface 
of  the  airfoil  does  not  involve  a complete  breakaway.  Particularly 
with  thin  airfoils  of  small  leading  edge  radius  the  separation  point 
is  "anchored"  at  the  leading  edge  followed  by  a reattachment  of  the 
separation- stream  line  at  a point  on  the  suction  surface  behind  the 
leading  edge.  A "bubble"  of  separated  flow  is  thus  formed  near  the 
leading  edge  within  which  the  velocities  are  quite  low  and  the 
perturbation  pressure  near  zero.  The  detailed  fluid  mechanical 
description  of  the  flow  is  very  complex,  but  the  net  effect  on  surface 
pressures  is  substancitally  as  described  above.  When  the  airfoils 
of  a cascade  are  configured  as  in  the  compressor  of  a gas  turbine  it 
is  thought  that  the  channelling  of  the  relative  flow  by  adjacent 
airfoils  inhibits  the  formation  of  separated  regions  extending  to 
the  trailing  edges  and  promotes  the  leading  edge  bubble  phenomenon. 

In  this  paper1  an  analytical  method  has  been  developed  to 

i 

predict  the  perturbed  aerodynamic  reactions  of  a flat  plate  airfoil 
with  leading  edge  bubble  of  variable  extent  undergoing  harmonic  os- 
cillatory motion.  This  technique  with  an  empirical  knowledge  of  the 
time  history  of  the  reattachment  point  can  be  used  to  predict  stall 
flutter.  Some  numerical  results  and  the  computer  plots  of  the  dynamic 
loops  are  presented. 
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Physical  Model 


A two  dimensional  thin  airfoil  with  small  camber  and  incidence 
is  considered.  The  flow  is  inviscid  and  incompressible.  The  prob- 
lem is  linear  due  to  (i)  small  time  dependent  displacements  (ii) 
small  or  zero  mean  angle  of  incidence  and  (iii)  small  perturbation 
velocities.  Under  these  conditions  a pressure  function  or  acceler- 
ation potential  p exists. 

As  a model,  a flat  airfoil  with  zero  thickness  on  the  x-axis 
is  taken.  In  Pig.  1 flow  is  parallel  to  the  x-axis.  The  flow 
separates  at  the  leading  edge  and  reattaches  at  a point  S on  the 
suction  surface.  The  time-dependent  positions  of  re-attachment  are 
shown  by  dotted  lines. 

Boundary  Conditions 

1.  The  perturbation  pressure  is  zero  in  the  bubble. 

y = 0”  0<x<2b 

y = 0+  S(t)<x<l 

4.  Complex  acceleration  potential  (W  = (Ji+iT)  is  continuous  at 
the  trailing  edge  (Kutta  condition)  and  the  separation  point. 

5.  Complex  acceleration  potential  vanishes  at  infinity  (i.e. 

W(Z)->  0 as  Z — >■  -<*>)  . 

6.  W (Z)  is  infinite  at  the  leading  edge.  (Integrable  singularity) 

7.  Harmonic  Oscillations.  _ 

(Displacements  are  given  by  h = he^wt  and  o=  ae-*wt) 

Mathematical  Development 

The  flow,  airfoil  and  the  boundary  conditions  are  transformed 
by  the  Conformal  transformation  C = £+  i'l  = . This  reduces  the 

difficulty  in  obtaining  the  solution  of  Laplace's  Equation.  Thus 
the  Poincare  boundary  value  problem  is  converted  to  the  solution  of 
a singular  integral  equation. 

The  prescribed  value  of  on  the  £-axis  is  given  by 

a n 

|*^  = 4b  ^ on  -1<  £<0  and  s(t)<f.<+l 

n=  0+ 


($=  0 in  the  cavity) 

2-  v-<lt + vfe»y 


3 a = ( — h v - ) v — 

-y  1 at  ax;  y ay 


where  s(t)=  ^S(t)/(2b) 

Fourier  transformation  is  used  to  find  the  solution  <P(£,nrt) 
of  the  Laplace  Equation  for  the  half  plane  problem  which  vanishes 
as  n— ► + w and  reduces  to  <|>(£»0+»t)  for  n=  0+. 

Using  Laplace's  Equation  and  applying  the  Fourier  transform 
wrt  the  variable  £ +00 


i#n»t)  =7*  l<K£,n»t)] 


e“3-*  d£ 


3 * (m,n*t). 


J 14> ] 


t* 

7 [<Dcc(£,n,t)]=  J< |,C5(£,n,t)e"jmC  d£ 


Using  partial  integration 

4,00 

7 -m2  J 4>(£,n,t)e~jmC  d£ 


a Q (mf  n , t)  2 


m 4>(m,n»t)=0 


The  solution  of  this  equation  is 
$(m,Ti*t)  =$  (m,0,t)e  ^ ^ 


Using  the  boundary  condition 

e-W"  [ Jut’.O't,.-™'  ari 


and  taking  the  inverse  Fourier  transform  of  4> 

4-eo 


♦ (£»n/t)  = 7_i[*  (m,n#  t)  ]■ 


*(m,n»t)e^m^  dm 


♦ (£#n,t)  - i /*♦(£»  0#t)  l f eljmU  ,m|1  dm]  d£ ' 


+ 0O 

U.n.t)  = i Jtu'.ott)  d5' 


This  is  Poisson's  integral  for.  a harmonic  function. 

The  conjugate  function  V is  found  using  the  Cauchy-Riemann  equation. 


-v(e,n,t) 


i [♦({lott)  ac‘-c(t) 

» ' n2+(E  -C) 


Using  the  indentity 


, r dy 

-YU.O?t)  = -¥<-l,0,t)  ~ J Tff' 


and  putting  n=  0+  in  the  previous  equation  results  in 


£ J<t>(z',ot 


t)  =» C(t) 

z'-z 


9 

= -Y(-l,o|t)  - J* 


The  range  of  integration  of  the  infinite  integral  may  be  reduced  as 
follows: 


i J <J>(£’/0tt)  * FU,t) 


where 


r 

F(?,t)  = - J <P(Z',0? t)  + w J 4>(Z’foU)  | 

-1  s(t) 

This  results  from  the  assumption  of  zero  perturbation  pressure  in  the 
separation  bubble  and  the  inability  of  a single  flat  airfoil  to  exhibit 
a pressure  perturbation  in  its  wake 

The  solution  of  this  last  integral  equation  for  the  two  sections 
between  (-1-rO)  and  (s(t)->  +1)  which  is  bounded  at  the  non-special 


HUM 


ends,  and  which  satisfies  the  singularity  condition  at  the  leading 
edge  and  Kutta  condition  at  the  trailing  is*4* 


where  the  constant  of  integration  C(t)  has  been  put  equal  to  zero 
since  -4* ->  0 as  n — > 00  • 


/ 


1 


where 


N{?)  - 


" " J tb  ar  " j *n 


Substituting  this  expression  in  Eq. (1)  one  obtains 


4>  (C,0+,  t)  = bX 


F (0  ) + Q(£) 


where  F 


C0)40«)--i 


~2  (-~  -'i'  (-1/  0+, t)  +N  (C  1 ) |^ 


J -»(-ifO\t)+N(C') 


s (t) 


Using  Eq. (2) 


-'P(-i,0+,t)+N(0)+Q(5)='i'(-i,0+,t)  ' | 

71  J l-r.2  5 * 


7 re 


-s  (t)  ) d£,  * 1 

,2  " e'-e  n 


1-e*  * ’ 


s(t) 
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1 _c  » ^ s,  “s 


s (t) 


Using  integration  over  the  contour  shown  in  Fig.  1 


e • <e • -s (t) ) dc  . | icu'-su))  d^v 

I — ^T2 — I J ,;;,2  r-c 


1-C 


s(t) 
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Using  the  identity  for  -'if  U , 0+,  t) , Eq.  (2)  and  the  expression  for 


N(£)  reveals  that 


FU,t)  = -YU,  0r,t) 


Assuming  harmonic  time  dependence,  the  Euler  equation 
3v  v 3v 3Y 


becomes 


av  + m 5 = _ i li 

3x  V V ax 


jut 

since  v=  v eJ 


Using  the  boundary  condition  v (-“)=  0,  the  above  differential 
equation  has  a solution  of  the  form 


_ 1 V 

v (x,0 ) = - ^ e V 


e V dx 
x 


where  the  path  of  integration  is  taken  to  lie  entirely  along  the 
real  axis.  The  above  equation  is  integrated  by  parts  to  reduce  the 


a 


order  of  the  singularity  in  the  integrand. 


v (0  + , 0 ) = - ^ 


- 

f xl 

- ? (0  + ,0)  + 

— V 

H<(x,0)  e 

_ 

dx 


(9) 


The  boundary  condition  is  applied  just  back  of  the  leading  edge. 


x=  0 


+e 


0 


Lira  | e 

e*  o 


j hOC 
V 


¥ (x, 0)dx=  Lim  j (1+ 
e -»-0 


^)  f?(0,n,t)ldx 

L n=+  /[? 


251 


+Lim 
e-y  0 


e 

J (i+J§5>  (x 


,0 


x_ 

2bl 


(10) 


Expressing  F(C/t)  as  a generalized  polynomial  in  £,  using  Eqs.(2), 
(5) , (6)  and  performing  the  integration,  it  is  found  that 


+ e 


Lim 

eVP 


3“  x 

e V V (x,  0)dx  = 0 


(ID 


-e 


Thus,  the  upper  limit  in  Eq. (9)  may  be  changed  to  0 . 
Writing  x=  — I x | 


v(0+,0)=  - ¥(0+,0)  + j VU, 


~3v  ,X| 

0)  e d)xl 


with 


o = 


= 13EL  and  ki? 


2b 


w 

Vv  (0+,  0)  = - *(0+.0)  + j2k  ^ 


- j2k 

e *(x,0)do 


Now 


, one  has  to  change  the  coordinates  of  the  function  ¥ . On  the 


lower  surface  x=0 ' , y-0  corresponds 


-8- 


r 


to  £ = 0“  , n = 0+  and  -»<x<o  > Y~Q 
corresponds  to  £=  0 , 0<n=  J frr  <“  . Thus,  continuing  the  previous 


reduction 


Vv(0  + ,0)=  - m,0  + ,t) 


+ j2k  J e~i2ka  (o,n,tj 


n=  o 


With  reference  to  Equations  (6)  and  (7),  changing  the  dummy  variable 
and  substituting  n=  a reduces  the  previous  equation  to 


co  0 

VV(0+,0)=F(0-)-3f*  /e-j2k°|  jjuf-s(t)) 


F ( 0, t)  +Q ( £)  da 


)2k  2kal  /j -i2 

TT  6 Y C (4-S  (t)) 

n s(t) 


F(0,t)+QU)]  ^ 

. ~ J o+t 


The  above  function  can  be  expressed  in  terms  of  tabulated 
functions 


s-j2ko 


— do  = g (2k£  ) - jf  (2kO 


where 


f (Z)  = sin  Z Ci (Z)  + cosZ(-Si(Z)) 

g(Z)  = sin  Z (-Si(Z))  - cosZ1  Ci(Z) 

3 5 7 

—Si  (Z)  = \ ~ z + fry  - - §7-5—  + frr  ' 


Ci  (Z)  = 0.577216 + JlnZ  - - + — 


2 ! 2 4 ! 4 6:6 


+ . . * 


Eq. (12)  is  solved  by  numerical  quadratures  yielding 


F(0,  t)  = 


Vv(0  + ,0)+j2k(R1+jJ1) 
l-j2k(RQ+jJo) 


-9 


where 


R +jJ„ 
o o 


1 

TT 


v. 


J 


luki 

\J  C-s 


S2) 


TtT 


(g-jf)dC  + 


i 


uirrju0- 


?-s(t) 


(g-jf)d? 


-1 


s (t) 


(14) 


and 


+1 


Rl+jJl= 


?< 


Q(S>  (g-:if>d5  + 


& 


Q(U  (g-jf)dQ 


-1 


s (t) 


Lift  and  Moment 

Taking  the  lift  to  be  positive  down; 
2b  2b 

L = - f p dx  - j (-p)dx 
0 s(t) 


where  the  conformal  transformation 

dx 


x = 2b(?2-n2) , 


d? 


= 4b? 


is  used  to  yield 


— 

0 

+1 

- 

Cl 

A 

<• 

1 

II 

14 

$(?»0)?d?  + 

r 

$(?,  0) ?d? 

Vj 

- 

- 

1 

s(t) 

Using 


r 0+/ t)—  j 


1-C 


[F  ( 0)  +Q  (?)  ] 


?<?-s(t)) 

which  was  derived  earlier  in  the  analysis. 


0 


+1 


L = -4bp 


[ f Si 


l F ( 0,  t)  +Q  (?,  t)l  d? 


* f iS 


(15) 


] 


^;s(|^-[F(0,t)+Q(?,t)]d? 


+1 


+ I ^ 


-s(t))(l-e*  ) 


V-Z 


+1' 


s (t) 


The  above  integrals  are  evaluated  by  contour  integration  (contours 
similar  to  Fig.  1).  When  the  result  is  substituted  in  the  previous 
equation  it  simplifies  to 


K* . si|U . s!|ti  + 1 , 


! 


**  JL  l-* 


(19) 


Fh(0>  = 2(  VjBh)b  » b 


(20) 


where  the  subscript  h denotes  bending  motion. 

_ 1 - 1 + 4k 

Ah+2Bh=  2 2 l-j2k(R0+jJ0)  (21) 


Ro+jJQ  is  given  by  Eq. (14) 

and 


s(t) 


Substituting  the  above  results  in  Eq.(l6),  lift  for  bending  motion 
is 


+ 


[ s (t)  € _ s ( t)  . -2  _j. 

1 — o n ^ 


2 2- 
2b  w h 


(23) 


Moment  for  bending  motion  is 


-8b2p 


^7^  } 2(Ah+jBh)b2w2h  + 2b2(o2h 


/T-s(t) 


l-Mt)_i„£l(tL+c2  + l]  ^ 


2(VjBh)bVh  + 2bV| 


s(t) 


r s(t)  -c 
2 


Torsional  Motion 


+ e2  ♦ b c2« 


For  torsion  about  the  leading  edge,  the  displacement  for 
a = a e^wt  positive  stalling,  is  expressed  as 

y = - axe^wt 

Analogously  to  the  bending  motion  case,  the  expressions 
are  given  successively 

v (C+) - -Va 


3$  r . 8 . 2 2 2 2r3,  - 

■^  = [-]  j b w 5 + 8b  u C 1 a 

N(<;)=  [-j  J bV  (C2-l)  + 2b2a>2U4-l)}  a 


Q(€f t)=  - - 


, .4,2  2,  - 
(~3j^b  u)  )a 


V (C,~s(t))(i-5,Z)  dC* 


T X 

f y5'(5' 

s(t) 


2 d£  • 

-s(t))(l-5’  ) frr^ 




-2bVs  f \v  w (c,-s(t))(i-c2)(iH,2)f^ 


C’-C 


J Jv{V- s(t))(l-5'2)  (l+c*2)|^  j.  -j  ^b2w2a+2b2aj25 

The  above  integrals  are  solved  by  similar  elaborate  contour 
integrations.  Substituting  the  results  and  after  some  algebra. 


r ...  2 2-  I .4  tr2  s (t 
Q(5,t)=b  2 


2 s (t)  g _s*(t)  11+- 

"2  8 2J  z 


/ i-4  s (t)  /I  s2(t),r2  .-s(t)  _ s3(t). , _ 5-s4_(tI  _ s?.(t)  3 t ^ 

j +{j  — 8 — n { 4 16  128  16.  8 

J (28) 

F(0-,t)  for  torsional  Motion  is  obtained  by  substituting  Eqs.(25),  (14), 
(15)  and  (28)  into  Eq.  (13) . 

Analogously  to  the  bending  motion,  the  lift  and  the  moment  for 
torsional  motion  are  obtained  by  substituting  the  above  values  of 
Q(£,t)  and  F (0 , t)  into  Eqs.  (16)  and  (17)  respectively. 

Special  Cases: 

(i)  Attached  Flow: 

For  Attached  Flow  s(t)  = 0. 

Using  this  value  in  Eq.  (1) 


*U,0+,t)=  - i 


,-,’2 


F(S’) 


-1 


The  above  result  satisfies  the  singularity  condition  at  the 
leading  edge  (5=0),  and  the  Kutta  condition  at  the  trailing 
edge(  +1). 

Eq. (19)  with  s(t)=  gives 

5u>-  2^5-  (?24> 

(ii)  Stalled  Airfoil: 


For  an  oscillating  airfoil  with  complete  separation  on  the  suction 
surface  s(t)  = 1. 

For  this  value  Eq. (1)  simplifies  to 

0 


({.o+,t)=  - i /If*  J r(C. 


t)  2il_ 
1 K'-K 


This  satisfies  the  singularity  and  Kutta  conditions  at  the 
leading  and  trailing  edges  respectively. 

For  full  separation,  Eq. (19)  of  bending  motion  becomes 

6(5)=,2bVh  (t2.1  ? + 3, 

and  Eq.  (28)  of  torsional  motion  simplifies  to 
z.,~\  ^ 2 , .2-  , ^4,^2  1 . 3 v ,n/r4  1 ,3,3^2  5 ^,35 

Q(C)=  b u>  o S+8>+2(*i  ~2  * +85  "T6  ^+l28  )] 


Similarly  the  lift  and  the  moment  for  these  special  cases  are 
obtained  by  substituting  the  relevant  expressions  derived  above  into 
Eqs.  (16)  and  (17) . 
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Stability  Analysis 


This  analysis  is  aimed  at  predicting  the  perturbed  aerodynamic 
characteristics  of  a flat  plate  airfoil  with  leading  edge  bubble 
separation.  The  motion  is  a harmonic  bending  and/or  torsional 
oscillation  in  an  otherwise  undisturbed  uniform  field.  The  flow 
separates  at  the  leading  edge  and  reattaches  at  a point  which  is 
moving  periodically  with  time  between  two  limits  on  the  suction 
surface. 

One  can  find  the  moment  coefficient  about  any  axis  using  the 
following  expression. 


- <Ie>  c, 


, x . 2 _ . x . „ 

«h0  Ma0  2b  \ ‘ 2b  Lh 


where  C = moment  coefficient  about  any  point  x along 

the  chord  due  to  torsional  oscillation  about 
this  point. 

CM  = moment  coefficient  about  the  leading  edge  due 
to  bending  oscillation. 

Cw  -~=  moment  coefficient  about  the  leading  edge  due 
M - 

a0  to  torsional  oscillation  about  the  leading  edge. 

CT  = lift  coefficient  due  to  bending  oscillation. 

C = lift  coefficient  due  to  torsional  oscillation 
a about  the  leading  edge. 

In  the  above  analysis,  the  point  of  reattachment  is  assumed  to 
be  given  by 


S (t) = ~ [Sj+Sj) +(S1-S2) cos (wt-x)l 


The  only  obstacle  to  applying  the  method  is  the  definition  of 
the  exact  time  dependency  of  the  point  of  reattachment.  The  present 

technique  is  semiempirical  in  the  sense  that’  it  uses  the  above 

/ 

expression  for  the  reattachment  point.  This  empiricism  in  specifying 
the  expression  can  be  improved  using  flow  visualization  techniques 
and  high  speed  photography  or  detailed  boundary  layer  analysis.  The 
variation  of  this  point  depends  on  a number  of  variables  including 


the  reduced  frequency  of  motion,  the  airfoil  geometry,  the  nature 
of  the  flow  pattern.  If  this  history  is  analyzed  for  different 
reduced  frequencies  from  experimental  data,  more  suitable  empirical 
functions  can  be  obtained. 

Given  this  time  dependent  movement  of  the  point,  the  present 
technique  can  indicate  the  occurrence  or  absence  of  stall  flutter  for 


any  configuration.  The  sign  of  the  value  of  work  done  per  cycle, 
positive,  zero  or  negative,  will  indicate  this  information  about 
stall  flutter. 

The  work  done  for  pure  bending  motion  is  represented  by 


where  CT  = CT  + C cos  (iot+4». ) + ... 

Lh  0 L1  1 

and  h = hQ  cos  tot 

Taking  the  integral 

W ^ - i it  CT  . sintj).. 

Some  numerical  values  are  presented  in  Table  2 as  examples  for 
the  coefficients  and  phase  angles  of  the  "Lift  Coefficient"  term 

Similarly  for  pure  torsional  motion 

2tr 

» ~ J c„  • at  • d<“t> 

0 

Numerical  Results 

A computer  program  has  been  developed  to  calculate  the  lift  and 
moment  coefficients  by  numerical  integration  using  the  expressions 
presented  above.  The  program  calculates  the  coefficients  at  a certain 
number  of  points  for  each  cycle  of  oscillation.  It  uses  harmonic 
analysis,  and  finds  the  coefficients  of  Fourier  series  which  rep- 
resents the  lift  and  the  moment  coefficients.  The  work  coefficient 


for  bending  oscillation  is  represented  by  one  of  these  Fourier 
constants , the  coefficient  of  the  first  "sine"  term.  The  Fourier 
series  is  then  employed  to  plot  the  graphs  of  the  lift  and  the 
moment  coefficient  loops.  The  graphs  presented  in  this  paper  are 
directly  obtained  using  this  computer  program  and  the  plotter  output. 

Results  and  Conclusions 

The  effect  on  the  lift  and  moment  of  a single  nonstationary 

airfoil  with  a leading  edge  bubble  of  variable  extent  described  in 

this  paper.  To  verify  that  the  developed  mathematical  model  is 

applicable,  fully  separated  (supercavitated)  flow,  which  is  a special 

case  of  the  general  method  presented  above,  is  studied  for  this 

purpose.  It  is  found  that  the  general  expressions  derived  for  lift 

and  moment  coefficients  simplify  exactly  to  the  results  obtained  by 

(2) 

previous  investigators  for  supercavitating  airfoils.  The  numerical 
results  are  also,  in  good  agreement  with  the  results  obtained  earlier 
as  shown  in  Table  1. 

An  initial  set  of  numerical  and  graphical  results  have  been 
obtained  for  pure  bending  (plunging)  oscillations  of  the  airfoil. 
Although  the  discussion  is  based  on  these  necessarily  preliminary 
data,  many  of  the  physical  conclusions  carry  over  to  expected  results 
to  be  obtained  in  the  future  for  pure  torsional  oscillations. 

The  location  of  the  bubble  that  forms  on  the  airfoil  governs 
the  airfoil  stall  characteristic.  Aerodynamic  reactions  vary  in  a 
periodic  but  non-sinusoidal  manner  for  an  airfoil  oscillating 
harmonically.  The  hysteresis  increases  as  the  travel  range  of  the 
reattachment  point  gets  wider.  It  is  believed  that  the  numerical 
results  deviate  somewhat  from  the  mathematically  exact  solution  if 
the  reattachment  point  approaches  the  immediate  neighborhood  of  the 


leading  edge  singularity.  The  amplitudes  of  lift  and  moment  due  to 
oscillation  become  smaller  as  the  point  of  reattachment  moves  closer 
to  the  trailing  edge.  This  is  attributable  to  the  reduction  of  the 
perturbation  pressure,  and  hence  the  contribution  to  unsteady  lift  and 
moment,  of  a larger  portion  of  the  suction  surface  including  the  leading 
edge. 

Damping  derivatives  (imaginary  part  of  lift,  CLh,  and  moment 
C„  , coefficients)  attain  larger  magnitudes  with  increasing  reduced 
frequencies.  These  higher  reduced  frequencies  produce  higher 
amplitudes  for  unsteady  aerodynamic  reactions  in  general,  see  Figs. 

6,7,8  and  Figs.  9 and  10. 

The  present  analysis  supplies  the  necesary  information  required 
to  determine  the  aerodynamic  work  including  its  sign.  The  work  can 
be  analytically  obtained  by  the  cyclic  integrals  ^Ldh  for  bending 
motion  and  ^Mda  for  torsional  motion.  From  the  type  of  data  shown 
in  Table  2 this  work  can  be  calculated  and  is  found  to  depend  on 
Cli  and  <{s^.  Aerodynamic  work  per  cycle  may  also  be  interpreted  as 
the  area  of  the  closed  figures  in  the  L,  h plane  or  the  M,  a plane. 

Work  done  per  cycle  depends  on  reduced  frequency  and  phase  lag  of  bubble 
movement  as  can  be  inferred  by  comparing  Figs.  3 to  6 for  phasing  effect 
and  Figs.  6 to  8 for  frequency  effect. 

Fig.  11  is  an  example  for  a fixed  point  of  reattachment  showing 
the  expected  elliptical  lift  loop.  The  hysteresis  effects  due  to 
bubble  movement  are  emphasized  in  Figs.  7 and  9 and  in  Figs.  8,10, 
and  12  in  which  the  reattachment  point  moves  during  the  cycle  of 
oscillation  for  each  series.  Crossing  of  the  lift  and  moment  loops 
has  been  reported  by  previous  investigators  as  a possible  effect  of 
stalling.  See  Figs.  3,  4,  5 and  13,  the  last  being  one  example  of  a 
moment  loop. 


It  can  be  concluded  that  the  integral  formulation  of  the 
mathematical  model  as  presented  here  is  a good  representation  of  the 
leading  edge  bubble  for  the  cases  discussed  above. 

With  a better  understanding  of  the  time  history  of  the 
reattachment  point  and  perhaps  using  a time  dependent  phase  lag, 
the  present  technique  can  be  studied  further  and  using  a -suitable 
mapping  function  can  be  applied  for  cascades  of  airfoils.  An 
important  next  step  is  the  execution  of  a numerical  and  graphical 
study  of  the  case  of  torsional  motion. 
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Table  1 

Fixed  Separation  Point  in  Translational  Oscillation 
(Fully  Separated  Case) 


Reduced  frequency 


Imaginary  Part  of 
Lift  Coefficient 


Present  L.C. Woods 


Imaginary  Part  of 
Moment  Coefficient 


.061 

.949 


.062 

.965 


Present 

.019 

.292 


L.C. Woods 

.019 

.302 


Table  2 

Translational  (Bending)  Oscillation 

h=h  co swt 
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CLo  + CLl  C°S  + cL2  cos  fc>t  +<t>2)  + 


fl 

fl 
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flo 
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-n/2 

-0.00581 

0.25164 

-88.3 
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84 

0.0323 

2/9 

1/9 
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Fig.  1 - Physical  model  and  corresponding  mathematical  model 
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Lift  coefficient  as  a function  of  dimensionless  bending 
velocity  at  x = n/3. 
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Lift  coefficient  as  a function  of  dimensionless  bending 
velocity  at  x = - it/2  and  k = 0.2 
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Fig.  10-  Lift  coefficient  as  a function  of  dimensionless  bending 
velocity  with  different  reattachment  point  movement  as 
compared  to  figure  8. 
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Fig.  13  - Example  of  moment  coefficient  loop  for  data  of  figure  8 
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